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We present a search for anomalous components of the quartic gauge boson coupling WWjj in 
events with an electron, a positron and missing transverse energy. The analyzed data correspond 
to 9.7 fb -1 of integrated luminosity collected by the DO detector in pp collisions at y/s = 1.96 TeV. 
The presence of anomalous quartic gauge couplings would manifest itself as an excess of boosted 
WW events. No such excess is found in the data, and we set the most stringent limits to date on 
the anomalous coupling parameters and . When a form factor with A cu tofi = 0.5 TeV is used, 
the observed upper limits at 95% C.L. are /A 2 \ < 0.0025 GeV -2 and \ag /A 2 | < 0.0092 GeV -2 . 
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I. INTRODUCTION 

In the standard model (SM) of particle physics, the 
couplings of fermions and gauge bosons are constrained 
by the gauge symmetries of the Lagrangian. The non- 
abelian gauge nature of the SM predicts the existence of 
trilinear (VVV) and quartic (FFW) gauge couplings 
(V = r y,W,Z). These include quartic couplings WW'yy 
between W bosons and photons that can be probed di- 
rectly at hadron colliders [IHIJ], but that are too small 
to be observed at the Tevatron, as will be shown later. 
Quartic couplingsprovide a window on electroweak sym- 
metry breaking [j, Q and can be probed by the mea- 
surement of W boson pair production via two photon 
exchange. 

Quartic couplings also allow for probing new physics 
that couples to electroweak bosons. As an example, the 
contribution of virtual heavy particles beyond the SM 
might manifest itself as a modification of the quartic cou- 
plings between W bosons and photons Observing 
the resulting anomalous couplings from such processes 
could be the first evidence of new physics in the elec- 
troweak sector of the SM. 

In this paper, we will focus on the search for WWjj 
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anomalous quartic gauge couplings (AQGCs) using data 
collected by the DO experiment at the Fermilab pp Teva- 
tron Collider, in events with an electron, a positron and 
missing transverse energy. The main production dia- 
grams are shown in Fig. [T] Pairs of W bosons are pro- 
duced via photon exchange, where the photons are di- 
rectly radiated from the colliding proton and antiproton. 
Triple gauge couplings WWj are assumed to be at their 
SM values (deviations from these values have been con- 
strained by the DO Collaboration Q and others (rol - [l3| ) . 

The parameterization of the AQGCs is based on 
Ref. [14(, and only the lowest dimension operators that 
have the correct Lorentz invariant structure and fulfill 
SU(2)c custodial symmetry [TH are considered. Such 
operators involving two W bosons and two photons are 
of dimension six: 

_ P 2 w 
~p 2 n w 

C$ = -^-^F m F^{W +a Wp +W~ a W+), (1) 

where F^ v is the electromagnetic field strength tensor 
and W^; is the W ± boson field, and are the usual 
notation for the parametrized quartic coupling constants, 
where a non-zero could be due to an exchange of a 
heavy neutral scalar, while heavy charged fermions would 
contribute to both off and Oq . The new scale A is in- 
troduced so that the Lagrangian density has the correct 
dimension of four and is interpreted as the typical mass 
scale of new physics. The current best 95% C.L. limits on 
these anomalous parameters come from the OPAL Col- 
laboration from measurement of WW"/, qqjj and vDjj 
production [l6[ at the CERN LEP Collider: 
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FIG. 1: Diagrams contributing to W boson pair production 
via photon exchange, with (a) triple WW"/ and (b) quartic 
VFW77 couplings. 



0.020 GeV~ 2 < a^/A 2 < 0.020 GeV" 
-0.052 GeV" 2 < a^/A 2 < 0.037 GeV~ 



(2) 



The pp — ¥ ppW + W~ cross section via photon ex- 
change rises quickly at high energies when the anomalous 
coupling parameters are non-zero, and manifests itself in 
particular with the production of boosted W boson pairs. 
In the SM, the 77 — > WW cross section is constant in the 
high-energy limit due to the cancellation between the rel- 
evant diagrams. When the new quartic terms are added, 
the cancellation does not hold and the cross section will 
grow to violate unitarity at high energies. This increase 
of the cross section can be regularized with a form factor 
that reduces the values of and Oq at high energy 
while not modifying them at lower energies. Following a 
standard approach, we introduce the following form fac- 
tor Q: 



.W 



,w 



(l + ^/Ac 2 utoff ) 2: 



(3) 



where M 77 is the invariant mass of the two photons, and 
A cu toff is chosen to be either 0.5 or 1 TeV, following the 
prescription of, e.g., Ref. [f|. In the following, we provide 
limits on anomalous couplings with and without form 
factors. 



II. DATA AND MONTE CARLO SAMPLES 

The full Run II set of data recorded by the DO detec- 
tor is considered in this analysis, representing 9.7 fb -1 of 
pp collisions at y/s — 1.96 TeV delivered by the Tevatron 
between 2002 and 2011, after the relevant data quality re- 
quirements are invoked. The DO detector used for Run II 
is described in detail in Ref. 



17j. The innermost part of 



the detector is composed of a central tracking system 
with a silicon microstrip tracker (SMT) and a central 
fiber tracker embedded within a 2 T solenoidal magnet. 



The tracking system is surrounded by a central preshower 
detector and a liquid- argon/uranium calorimeter with 
electromagnetic, fine, and coarse hadronic sections. The 
central calorimeter (CC) covers pseudorapidity [lH \Vd\ 
< 1.1. Two end calorimeters (EC) extend the coverage to 
1-4 ^ \Wd\ ^ 4.2. Energy sampling in the region between 
the ECs and CC is improved by the addition of scintil- 
lating tiles. A muon spectrometer, with pseudorapidity 
coverage of \r)d\ < 2, resides outside the calorimetry and 
is comprised of drift tubes, scintillation counters, and 
toroidal magnets. Trigger decisions are based on informa- 
tion from the tracking detectors, calorimeters, and muon 
spectrometer. Details on the reconstruction and identifi- 
cation criteria for electrons, jets, and missing transverse 
energy, ]fir, can be found elsewhere [l9| . In this paper 
we call both electrons and positrons "electrons," with the 
charge of the particle determined from the curvature of 
the associated tracks in the central tracking system. 

The background where, like the signal, the proton and 
the antiproton are intact in the final state, originates 
from photon exchange and double pomeron exchange 
(DPE) processes [2(|. Both these backgrounds and the 
AQGC signals are modeled using the FPMC [2l[ genera- 
tor, followed by a detailed GEANT3-based [22j simulation 
of the DO detector. Data from random beam crossings 
are overlaid on the MC events to account for detector 
noise and additional pp interactions. The predictions of 
the FPMC generator, which are made assuming that the 
proton and antiproton are left intact after the interac- 
tion, are consistent with those of the lpair [23| genera- 
tor, which in turn are consistent with the measurement 
of the cross section for exclusive e + e~ production by the 
CDF Collaboration [Hj]. 

Diffractive and photon exchange backgrounds to this 
search are exclusive e + e~ and t + t~ production through 
i-channel photon exchange (Drell-Yan) and inclusive 
W + W~, e + e~, and t + t~ production through DPE. 

Since the outgoing intact proton and antiproton are 
not detected in this measurement, we also need to con- 
sider non-diffractive backgrounds. These backgrounds 
are Z/ 7 *+jets, tt and diboson (W + W~ , W ± Z and ZZ) 
production, and processes in which jets are misidentified 
as electrons: FF+jets and multijet production. The sim- 
ulated samples used to model them are identical to those 
described in Ref. [l!|. All of these backgrounds, except 
multijet production, are modeled using the PYTHIA [25| 
or ALPGEN [26j generator, with pythia providing show- 
ering and hadronization in the latter case, using the 
CTEQ6L1 H3 parton distribution functions (PDFs). 
The multijet background is determined from the data by 
inver ting some electron selection criteria, as described in 
Ref. [lj. 

Single diffractive (SD) processes, for which either the 
incoming proton or antiproton is intact after the inter- 
action while the other is destroyed, have similar features 
to non-diffractive (ND) processes in the direction of the 
broken proton or antiproton, contrary to DPE processes 
where both the proton and antiproton are intact. Since 
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the cross section ratio of SD to ND processes is about 
(2-3)%, which is below the uncertainty on cross sections 
of ND processes cross sections, the contribution of SD 
processes is neglected in this analysis. 

The selection of data events is similar but more strict 
than the search for the Higgs boson in the H — » 
W + W~ — > e + ve~v channel that is described in de- 
tail elsewhere which includes the same trigger ap- 
proach with no explicit requirement. A preselection is 
applied to the data by requiring two high-transverse mo- 
mentum (high-pr) electrons with opposite charge. The 
leading- and trailing-px electrons are required to satisfy 
pf^ > 15GeV and p^ > lOGeV, and their invariant 
mass is required to be M ee > 15 GeV. In addition, these 
electrons are required to be within the acceptance of the 
calorimeter (|%| < 1.1 or 1.5 < \r) d \ < 2.5 [18]), with 
at least one electron required to be in the central part of 
the calorimeter (1^1 < 1-1). The only difference from the 
event selection in the Higgs boson search is that we veto 
events with at least one jet with pt > 20 GeV, \r/d\ < 2.4, 
and matched to at least two tracks associated with the pp 
interaction vertex. The inclusive cross section for exclu- 
sive W boson pair production through photon exchange 
in the SM at Vs = 2 TeV is cr(pp -> ppWW) = 3 fb, but 
after the preselection only 0.1 event is expected from this 
process, unless it is enhanced by AQGCs. 

To correct for any possible mismodeling of the lepton 
reconstruction and trigger efficiencies, and to reduce the 
impact of the luminosity uncertainty, scale factors are 
applied to the Monte Carlo (MC) samples at the prese- 
lection stage to match the data. The Z boson mass peak 
region in the data and MC samples after the preselec- 
tion is used to determine normalization factors. Their 
differences from unity are found to be consistent with 
the luminosity uncertainty of 6.1% [281 ] . The pr distri- 
bution of Z bosons is weighted to match the distribution 
observed in data [29| , taking into account its dependence 
on the number of reconstructed jets. The pt distribution 
of W bosons is weighted to match the measured Z boson 
Pt spectrum, corrected for the differences between the 
W and Z p T spectra predicted in NNLO QCD [3(j. The 
distribution of the pr of the leading electron after the 
preselection is shown in Fig. [Ha). 

Following the same strategy as described in Ref. fljjj . 
boosted decision trees (BDT) are used to reject the large 
Z/j*+jets background, that is dominant after the pre- 
selection. The input variables to this "selection BDT" 
are kinematic quantities, including the electron mo- 
menta, the azimuthal opening angle between the two 
electrons, $t, variables that take into account both 
]fir and its direction relative to each electron, and ob- 
servables that differentiate between real and misrecon- 
structed Ipr- The cut on the selection BDT, which de- 
fines the final selection, is chosen such that the contri- 
butions of the Z/7*+jcts, W+jets, and diboson back- 
grounds are of comparable magnitude. The distribution 
of the single most discriminating variable, the transverse 
mass of the Ipr and the dielectron pair (Mr(ee, Ipr) = 



y/2 ■ p?^ ■ $t ■ [1 — cos A<j)(ee, $r)]), after the final selec- 
tion is shown in Fig. &[b). The expected and observed 
numbers of events after the preselection and the final se- 
lection are given in Table HI 

A final BDT is trained to separate the AQGC signal 
from all the other backgrounds. The same BDT is used in 
the study of both parameters and a^, , which feature 
identical kinematic characteristics. This BDT relies on 
the input variables of the selection BDT, complemented 
with additional variables characterizing the electron re- 
construction quality to discriminate against the instru- 
mental backgrounds (multijet and VF+jets production). 
The distribution of the final BDT output is shown in 
Fig. [2fc) and demonstrates the good agreement between 
the data and the background expectation. 

III. SYSTEMATIC UNCERTAINTIES 

Systematic uncertainties are estimated for the signal 
and for each background process. They can affect only 
the normalization or both the normalization and the 
shape of the final discriminant. 

Sources of systematic uncertainty that affect only the 
normalization arise from the uncertainties on the the- 
oretical cross sections of Z+jets (6%), W+jets (16%), 
diboson (6%), and ti (7%) processes; the multijet nor- 
malization (30%); and the modeling of the Iftx for the 
Z+jets background (5%). The diffractive backgrounds 
have been assigned a 100% uncertainty on their cross sec- 
tions due to the large uncertainties on the gluon density 
(for processes induced by pomeron exchange; the uncer- 
tainty on the gluon density inside the pomeron can reach 
40%, translating into an uncertainty of a factor up to 
2 on the cross section) and on the proton dissociation 
(for processes induced by photon exchange). For the lat- 
ter process, a 20% uncertainty has been assigned to the 
signal theoretical cross section. 

The sources of systematic uncertainty that also affect 
the shape of the final discriminant distribution are quoted 
here as average fractional uncertainty across bins of the 
final discriminant distribution for all backgrounds: jet 
energy scale (4%), jet resolution (0.5%), Ifrr modeling 
(4%), jet identification (2%), jet association to the hard- 
scatter primary pp interaction vertex (2%), and H^+jets 
modeling (10%). The systematic uncertainties due to 
the modeling of the px{WW) and the Acf> between the 
leptons, and the pr of the vector boson from the W+jets 
and Z+jets production (see Ref. [l!|) are less than 1% 
and taken into account. 



IV. RESULTS 

The data are found to be in good agreement with the 
background-only prediction, and upper limits are set on 
the anomalous parameters and . The modified fre- 
quentist CL S method (3lj is employed to set limits on the 
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P? 1 [GeV] M T (ee,£ T ) [GeV] Final Discriminant 



FIG. 2: [color online] The (a) leading electron pr at the preselection level, (b) the transverse mass of the Jfir and the two 
electrons after the final selection, and (c) the output of the final BDT discriminant after the final selection. In (a) and (b), the 
last bin includes all events above the upper bound of the histogram. The hatched bands show the total systematic uncertainty 
on the background prediction, and the signal distributions are those expected for af / K 2 = 5 x 10~ 4 GeV~ 2 and no form factor. 



TABLE I: Observed and expected numbers of events after the preselection and the final selection for data, signal (aJ^/A 2 = 
5 x ICP 4 GeV -2 and no form factor), and the different backgrounds considered in the analysis ("Diff" stands for the diffractive 
backgrounds). 





Data 


Total background 


Signal 


Z/"* — > ee 


Z/'y* — > tt 


tt 


VF+jets 


Diboson 


Multijet 


Diff. 


Preselection: 


572700 


576576 ± 11532 


12.2 


566800 


4726 


15 


623 


517 


2716 


1180 


Final selection: 


946 


983 ± 108 


11.6 


291 


22 


8 


370 


287 


5.4 


0.2 



AQGCs, where the test statistic is a log-likelihood ratio 
(LLR) for the background-only and signal+background 
hypotheses. The LLR is obtained by summing the LLR 
values of the bins of the final BDT output. In the 
LLR calculation, the signal and background rates are 
functions of the systematic uncertainties that are taken 
into account as nuisance parameters with Gaussian pri- 
ors. Their degrading effect is reduced by fitting the 
expected contributions to the data by maximizing the 
profile likelihood function for the background-only and 
signal+background hypotheses separately, appropriately 
taking into account all correlations between the system- 
atic uncertainties [32j |. 

The 95% C.L. allowed ranges for the anomalous param- 
eter {etc) can be found in Table HT1 (IIIIll . assuming a^, 
(a,Q ) is zero. The limits are quoted both without a form 
factor and for a form factor with A cuto ff = 1 or 0.5 TeV (as 
advised, e.g., in Ref. @). The two-parameter limits are 
shown in Fig. [3] for different assumptions about the sig- 
nal, namely if no form factor is used and if a form factor 
is used with A cuto ff = 1 or 0.5 TeV. The two-parameter 
68% C.L. (95% C.L.) limits define the range of values 
of the anomalous coupling parameters for which the the- 
oretical cross section is lower than the upper 68% C.L. 
(95% C.L.) limit on the signal cross section, obtained in 
the single parameter limits. The effect of the presence 
of a Higgs boson with Mh = 125 GeV is not accounted 
for, but is expected to contribute less than 4 events af- 
ter the final selection, having kinematic distributions dis- 



tinct from signal, and to broaden the allowed ranges for 
the anomalous parameters by a negligible amount. 



V. CONCLUSION 

We have searched for anomalous WWj~f quartic gauge 
boson couplings by analyzing 9.7 fb _1 of integrated lu- 
minosity in the W + W~ — > e + ve~v final state using the 
DO detector. No excess above the background expecta- 
tion has been found. When a form factor with A cuto ff = 
0.5 TeV is used, the observed upper limits at 95% C.L. are 
\a^/A 2 \ < 0.0025 GeV" 2 and \a™/A 2 \ < 0.0092 GeV" 2 . 
These are a factor 4 to 8 more stringent constraints on 
and than the previous limits 16], and the only 
published limits to date from a hadron collider. 
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stitutions, and acknowledge support from the DOE and 
NSF (USA); CEA and CNRS/IN2P3 (France); MON, 
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ciencias (Colombia); CONACyT (Mexico); NRF (Ko- 
rea); FOM (The Netherlands); STFC and the Royal So- 
ciety (United Kingdom); MSMT and GACR (Czech Re- 
public); BMBF and DFG (Germany); SFI (Ireland); The 
Swedish Research Council (Sweden); and CAS and CNSF 
(China). 
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TABLE II: Expected and observed 95% C.L upper limits on |a,^/A 2 j 
the form factor. 


assuming a™ 


is zero and for different assumptions about 


Cutoff Expected upper limit [GcV~ 




Observed upper limit [GeV -2 ] 


No form factor 0.00043 
A cutoff = 1 TeV 0.00092 
A cutoft = 0.5 TeV 0.0025 




0.00043 
0.00089 
0.0025 


TABLE III: Expected and observed 95% C.L upper limits on \oc /A 2 | 
the form factor. 


, assuming Oq 


is zero and for different assumptions about 


Cutoff Expected upper limit [GcV~ 




Observed upper limit [GeV -2 ] 


No form factor 0.0016 
A cutoB = 1 TeV 0.0033 
A cut off = 0.5 TeV 0.0090 




0.0015 
0.0033 
0.0092 
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95% C.L. Contour 
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0) 
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CM 
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M " 0.03 r 

> 
at 

— 0.02 

CM 

< 
So 
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DO, 9.7 fb\ ee + P T 
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Data 
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95% C.L. Contour 



-0.01 
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-0.02 



0.005 
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FIG. 3: [color online] Two-parameter 68% and 95% C.L limits with different assumptions about the signal: (a) no form factor, 
or a form factor with (b) A cuto fi = 1 or (c) 0.5 TeV. 
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